Undulation instability of lipid membranes under an electric field 
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Abstract 

The influence of an electric field on a poorly conductive membrane such as a 
lipid bilayer is studied theoretically. The unbalanced electric stress created by 
an ionic current across a non-perfectly flat membrane gives rise to a destabi- 
lizing surface energy enhancing undulations. The deformation of a membrane 
attached to a frame and the subsequent force on the frame are derived and 
the electrohydrodynamic instability of a free floating membrane is also stud- 
ied. We find a most unstable mode of undulation, of wavelength in the [im 
range, connected to the crossover between membrane and solvent dominated 
dissipations. 
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Due to their low permeability to electrolytes, biological or other lipid membranes are 
strongly influenced by applied electric fields. Researchers are actively investigating phe- 
nomena such as electroporation (creation of long-lived pores in a lipid membrane under the 
action of a strong electric field) , in part because of very promising applications for gene 
therapy. Similarly, the electro formation of lipid vesicles is a widely used technique to form 
large unilamellar vesicles under the action of an electric field Despite its importance, 
there is no clear understanding of this experimental tool, which remains mostly empirical. 
Many other situations (electrofusion electro-injection of macromolecules in vesicles and 
cells ||) require a better understanding of the action of electric fields on lipid bilayers. This 
paper addresses this issue theoretically, and discusses the effect of an electric field on i) the 
force exerted on a rigid frame by an attached membrane, and ii) the dynamical instability 
of a free floating membrane and its relevance for the electroformation of neutral liposomes. 
Before reaching these two topics, we study the stability of an almost flat membrane under 
electric field. 

The force acting on an interface in an electric field E can be calculated by evaluating the 
discontinuity of the Maxwell stress tensor = e(EiEj — jE 2 5ij) across the interface 0. 
We consider an infinite membrane of dielectric constant e m (~ 2eo) and conductivity Xm"(— 
10~ 6 S/m) in a solvent of dielectric constant ex(— 80e ) and conductivity Xi {— lO^S/m) 
(numbers are typical for a lipid membrane in water, eo is the permitivity of vacuum). A 
fixed electric field E (in practice, a fixed electric current j) is applied across the membrane, 
and we note E m the electric field inside the membrane. 

If the membrane is perfectly flat, the electric stress is symmetrically balanced on both 
sides of the membrane (i.e., a e i ~ e m £^-6 m £^ = 0), which experiences a mere compression 
0. As already pointed out in 0, this perfect cancellation of the net electric stress is 
however accidental since any membrane curvature leads to an unbalanced net stress whose 
order of magnitude can be expressed as a e i ~ emE^d/R, where d is the thickness of the 
membrane and R its local radius of curvature (a curvature of the membrane leads to a 
non- uniform surface charge density - see Fig.l). In the case of a closed vesicle @, E m is 
given by the potential drop at the vesicle scale (i.e., E m ~ ER/d) as most of the electric 
field goes around the finite sized object, whereas the current continuity across the (infinite) 
membrane dictates the expression of E m here, namely E m ~ Exm/Xi- This resulting stress 
o e \ tends to enhance the (small) local undulations of the membrane u(r) = ^2 q u q e iqr (Monge 
representation in real and Fourier space). After a short transient time r e ~ (e m /d)xiH ~ 
10 -6 — 10 _5 s (H is the distance between the electrodes), the solution of the Poisson equation 
(in the limit e m Xm ^ e iXi) yields the following expression for a given Fourier mode: a e i = 
2e m (Exm/ Xi) 2 ( eqd ~ l)/(e qd + l)tt(r). Hence, integrating the work of the electric stress (in 
the large wavelength limit qd <C 1) yields a net decrease in energy 



This amounts to an effective negative surface tension T e i acting on the membrane. In the 
presence of a lipid reservoir (or at the interface between two immiscible fluids of different 
conductivities), this term leads to an electric field induced decrease of the interfacial tension, 
known as the electrocapillary effect ||. On the contrary, for a fixed number of surfactants, 
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enhancing membrane undulations under electric field builds a concomitant mechanical ten- 
sion because of deviation from the nominal area per molecule. Hence we expect a tense 
yet floppy-looking membrane in this case. For typical values of the electric field used in 
the electroporation and electroformation experiments E ~ 10 3 V/m, the electrostatic surface 
tension reaches T e i = 1CT 3 J/m 2 , which is of order the mechanical tension needed to rupture 
a lipid membrane [||]. 

Two consequences of this destabilizing effect are studied below: i) a membrane on a fixed 
frame is deformed until the electric stress is balanced by an opposing mechanical stress and 
ii) a free membrane undergoes strong undulations under electric field. 

o Case of a bilayer attached to a fixed frame 

The elastic behavior of a lipid membrane is generally characterized by a bending modulus 
k(~5x 10~ 20 J), and a stretching modulus K s (~ 0.1 J/m 2 ) |J. Typically, bending a lipid 
membrane involves energies of order 1 — 10k bT, while stretching it requires much larger 
energies. To simplify the description below, we will not include the bending rigidity in 
the treatment of the static deformation of the membrane, nor will we include the thermal 
fluctuations of the membrane, restricting ourselves to large electric fields. 

We write a linear theory for small deformations of the membrane qu q <C 1. The mem- 
brane is characterized by its total area S, its projected area S p , and its optimum area for 
which it is not stretched So- We use the two small parameters 5 = (S p — So) /So and 
X = (S — Sp)/So — (1 + 5)\ Yl Q 2 \ u q\ 2 — \ Yl Q 2 \ u q\ 2 - The electric-field-induced undulation 
term involves the parameter X only (Eq. ( P-l|) ), and the stretching energy involves the total 
increase of area with respect to So'. 

F = S \^K S {5 + 1) 2 - r e ,j) (0.2) 

Note that the electric- field-induced undulation term (T e i) is quadratic in the small quantity 
qu q , while the stretching term (K s ) corresponds to an expansion to the fourth order, con- 
sistent with the large value of the ratio of the stretching over the electrostatic parameters 
(3 = K s /Y el > 100. 

The minimization of this energy with respect to the out-of-plane membrane deformation 
X leads to an equilibrium class of membrane shapes corresponding to a total area difference 
A s 

A 5 = 5 + X e9 = ^~10- 2 (0.3) 
i r 2 

The energy of the membrane only is F m = ^Sq-^ 1 , and the energy of the system, including 

the electrostatic energy, is F = F m + F ei = — ^Sq-^ 1 + S T e i8. 

The surface tension of the lipid membrane at equilibrium is given by 7 = dF m /dS with 
S = S (l + 6 + X) 

At equilibrium, the mechanical tension is equal, but opposite in sign, to the electrostatic 
energy per unit area: 7 = T e i. The building of a mechanical tension in a membrane under 
electric field is due to an increase of its area. 
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Beyond the electroporation threshold, holes in the membrane reduce the mechanical 
tension by decreasing the membrane area and strongly increasing the membrane electrical 
conductivity 1/xm- The effect of the electric field could be measured via the setup depicted 
in Fig.2. The force measured by the spring is / = dF tot /dL p = T e iL y (~ iO~ 3 iV/m), where 
L p and L y are the membrane projected size in the direction of and perpendicular to the 
spring respectively. The membrane pulls on the spring (/ > 0) if I > —5. The force in the 
absence of field is simply: / = K s L y 5. 

It is clear from Eq. ( |0.3| ) that the balance between electric-field-induced undulation and 
membrane stretching does not select a particular equilibrium membrane shape, as both 
effects depend on the global increase of membrane area only. Including the bending energy 
of the membrane selects the shape of lowest curvature, namely the first harmonic q\ = tt/L p . 
The energy gap between the different harmonics is however small (of order ksT) for qX K < 1, 
where X K = 27Ta/ n/T e i. For large electric fields (E ~ 10 3 V/m), X K is very small (~ lOnm). 
For small electric fields, the thermal fluctuations dominate both the membrane shape and 
tension. We postpone the study of the interesting crossover between these two limits to a 
future publication. 

There are some dynamical issues connected to force measurements on lipid membranes 
(Fig.2). Lipid membranes spread on a frame are generally connected to the frame by meniscii 
which are large compared to the membrane thickness, and from which lipid molecules may 
flow toward the deformed membrane (a phenomenon reminiscent of the Marangoni effect 
|1(J). In this case, the membrane elastic stress consequent to the action of an electric field is 



rapidly released by migration of lipid molecules, and might not be detectable experimentally. 
The kind of force measurement depicted on Fig.l should however be feasible on polymeric 



membranes, which have shown sensitivity to electric field as well |lT[, and for which the 
migration of molecules along the membrane, if any, is very slow. An alternative experimental 
setup for lipid membranes involves measurements on lipid membrane strongly bound to a 
solid substrate, which are expected to exhibit a fairly slow lateral dynamics. 

o Instability of free floating membranes under E-field 

We now study the case of a free floating membrane under an electric field. This study is 
partly motivated by the experimental technique of electroformation of vesicles, which allows 
a controlled swelling of an electrode deposit of lipids to form vesicles of fairly well controlled 
sizes under electric field ||. Interestingly, this technique works for charged and nonionic 
lipids, although the optimal conditions vary widely with the nature of the lipids. 

The dynamics of a free membrane under electric field can be decomposed into two mech- 
anisms, i) The normal deformation of the membrane (described in the previous section), 
which creates a tension 7 in the membrane and saturates when 7 = T e i. ii) The lateral slid- 
ing (contraction) of the membrane, in an attempt to release this tension. The characteristic 
time for the normal deformation for a given undulation mode of wavelength A is rj_ ~ i]X/T e i 
(~ lCT 6 s for A = lfim), while the sliding motion involves the whole membrane t\\ ~ i]L/T e i 
(~ 10 _3 s for L ~ 1mm). The two mechanisms occur at very different timescales, and can 
be treated separately. 

A thorough treatment of the undulation modes of a film immersed in a solvent can be 



found in the literature [12]. The electrohydrodynamic instability of a layer of non-conducting 



fluid between two semi-infinite conducting fluids has been studied in [|T|], where special 
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attention is given to the peristaltic deformation modes (the two interfaces undulating in 
antiphase), as these modes lead to the destruction of the film when the two interfaces make 
contact. In the case of a lipid membrane, these peristaltic modes are suppressed because of 
the very low compressibility of the film. We study below the bending instability (interfaces 
undulating in phase) of the membrane. 

i) The normal displacement of the membrane involves viscous dissipation in and around 
the membrane jnyr5J. There are three main dissipation mechanisms, namely the dissipation 
in the solvent (viscosity rj = 10~ 3 Pa. s), which dominates the dynamics of large wavelengths 
deformation, the friction between the two monolayers (friction coefficient bf r — 10 8 Pa.s/m), 
dominant at intermediate wavelengths, and the membrane surface dissipation (surface vis- 
cosity fi = 10~ 10 P 'a. s.m). For undulating membranes, the latter mechanism is relevant 
at very small wavelength of the order of the bilayer thickness d = 5nm only, and will be 
neglected in what follows. We present below a simplified description of the interplay be- 
tween external and internal dynamics. For a thorough treatment of membrane dynamics, 



see Seifert [14 



Neglecting inertia, a normal deformation of the membrane (of typical lateral size 2ir/q 
and typical velocity u q ) creates a motion in the surrounding fluid which propagates to a 
distance ~ 1/q. The curvature of the membrane leads to a velocity difference of order 
Sv = qdu q between the two monolayers. The power dissipated by viscous effect around and 
in the membrane can be written respectively: 

P v = v I dV (Vv) 2 = SriJ2 ?«J 

P fr = b fr ! dS {8vf = Sbf r A 2 (0-4) 
J 1 

This dissipated power must compensate the power stored in the membrane P m = dtF (the 
energy F is given by Eq.( p.2| )). This condition leads to an evolution equation for each 
deformation mode: 

{<qq + b fr d 2 q 2 )u q {t) = 

[T el q 2 (1 - f3(5(t) + At))) - ^] u q (t) (0.5) 

with (3 = K s /T e i > 100. The left hand side of this equation describes the viscous dissipation 
in and around the membrane and the right hand side consists of the electric-field-induced 
undulation term (r e /) with a stretching saturation (/3 term). Note that the bending rigidity 
of the membrane has been added to the membrane energy (k term), for it is mandatory for 
the description of the small wavelength deformations q > y/Y e i/ k. This equation is non- 
linear, since the saturation involves the total increase of area I = | 'Yl, q q 2 \u q \ 2 . Thanks to 
the different time scales for normal and lateral motions, the projected area of the membrane 
Spit) oc 5 (t) can be considered as constant for the short time evolution. 

This non-linear equation can be solved, but for our purpose, it is sufficient to consider its 
linearized form, and to treat the saturation dynamics separately. The short time evolution 
is described by the linear equation 

(rjq + b fr d 2 q 2 )u q = (T e iq 2 - nq A )u q (0.6) 
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The amplitude of a given Fourier mode has a time evolution u q (t) ~ e aqt with 



a i = —<l\ 1 , TiZ (°- 7 ) 



v V 1 + q/qdyn 

with the two characteristic wavevectors: 

q s tat = \f^r~5.l0 7 m- 1 qdyn = 1 !L~5.l0 5 m- 1 (0.8) 
The evolution rate presents a sharp maximum at 

? t = W w ) 1/3 = 8.10 6 ^ 1 (0-9) 

for q stat ^> g dj ,„. This defines a particular lengthscale which grows exponentially faster than 
the others: A* = 2n/q* ~ 1 /im. The corresponding growth rate is a* ~ ~ S.IO 5 ^" 1 . The 
evolution saturates similarly for all lengthscales when /3 [5(t) + Z(t)] = 1, at which point a 
mechanical tension T e i is established, leading to a contraction of the membrane. 

ii) The lateral motion of the membrane occurs at an almost constant surface tension, be- 
cause the time needed to build up the tension t± is much smaller than the time t\\ over which 
it could be released. Since the contraction of the membrane involves solvent flow over large 



lengthscales (mm), inertial effects must be included ||16|| . The flow created by the lateral mo- 
tion of the membrane (of size L) propagates over a size L z = L/(a/1 + L 2 / {vt)) [y = r//p is 
the kinematic viscosity ~ 10 _6 m 2 /s for water - p is the density of water). The comparison of 
the power dissipated by the sliding motion of velocity L: Pdi SS = V f dV{S/v) 2 = r]L 2 L 2 /L z , 
with the power stored in the membrane P s tor — ^eiLL, leads to the dynamical equation for 
the membrane size: 

Tp/ L r ITp/ 1 . 

L= ~-^ = ; = 0.10 

v l rj V /TTX7K) 

The short time evolution (t < L 2 =0 / v ~ Is) is dominated by the diffusion of the solvent 
flow. Assuming that L t =o = L Q (No stretching without electric field), the evolution equation 
for short time is 



L(t) = L ^l-A/3^ty(r 2 r u ) 
r r = vLo/Tei ^ 10~ 3 s r v = L 2 Q /u ~ Is (0.11) 

from which emerges a characteristic time scale T s ud e = ( T r T iy) 1 ^ 3 ~ 10 _3 s. This result is 
reminiscent of the dynamics of bursting of a soap film in a viscous environment JlT]]. The 
evolution of the variable 5 = L p /L — 1 which appears in Eq. (p.5|) , is given by 5(t) ~ 
2/3(t/r sMe ) 3/2 . 

In many practical situations, fixed boundaries or walls (the electrodes in the case of 
the electroformation of vesicles) may strongly modified the hydrodynamics. The solvent 
viscous dissipation (Eq. (|0~4]) ) is modified by the presence of a wall at a distance h from the 
membrane, and becomes P^h = Sr]J2 q u 2 /(q 2 h 3 ) for qh < 1. The left hand side of Eq. (p.6[) 
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is modified accordingly, and shows an optimal wavevector ~ \Jq* /h where q* is given 
by Eq.( |0.9| ). For h = lOnm, the fastest growing wavelength is of order 2k /q^ ~ 0.1/wn. 
The sliding motion is strongly affected by the wall. For times larger than h 2 /u, the time 
dependent length L z should be replaced by the distance h in Eq.( O.lOj ), and the dynamics 
is described by L(t) = L a/1 — tjrh with = TrL /h ~ 10 2 s instead of T s ud e = 10~ 3 s (note 
that solvent permeation through lipid membrane becomes relevant near a wall)! 

To conclude this paper, we would like to propose that the electrohydrodynamic instability 
described above may play an important role in the first stage of the electroformation of 
liposomes. A remarkable feature of this technique is that it produces vesicles of fairly 
well defined size. We find a fastest growing undulation mode of wavelength in the fim 
range (Eq. ( p.8| , p.9| )). This mode might be the precurssor of large scale deformations of the 
membrane which, after a complex process partly sketched in Ref || and involving coalescence 
of neighboring blisters, lead to the formation of closed vesicles (the size of which can reach 
50 fim for nonionic lipids |1§)- Future developments will include the treatment of small 
pores which are expected to be present in a membrane under tension, and their influence on 
both the membrane electrical conductivity (hence r e ;) and dynamics (solvent permeation 
through the membrane). 

We would like to thank A. Johner, M. Angelova, R. Bruinsma, J.F. Joanny, A. Ajdari 
and J.B Fournier for stimulating discussion and useful comments. 
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Fig.l Net accumulation of conduction charges near 
a curved lipid membrane under electric field 




Fig. 2. Possible experimental setup for the force measurement 
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